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We discuss how to extract information about the cosmological con- 
stant from the Wheeler-DeWitt equation, considered as an eigenvalue of a 
Sturm-Liouville problem. The equation is approximated to one loop with 
the help of a variational approach with Gaussian trial wave functionals. A 
canonical decomposition of modes is used to separate transverse-traceless 
tensors (graviton) from ghosts and scalar. We show that no ghosts appear 
in the final evaluation of the cosmological constant. A zeta function reg- 
ularization is used to handle with divergences. A renormalization proce- 
dure is introduced to remove the infinities together with a renormalization 
group equation. A brief discussion on the extension to a f(R) theory is 
considered. 

1 Introduction 

The Friedmann- Robertson- Walker model of the universe, based on the Einstein's 
field equations gives an explanation of why the Universe is in an acceleration 
phase. This is supported by data observations on type I supernovae[T]. Never- 
theless, to obtain such an expansion we need almost 76% of what is known as 
Dark Energy. Dark Energy is based on the following equation of state P = ujp 
(where P and p are the pressure of the fluid and the energy density, respec- 
tively). When lu < —1/3, we are in the Dark energy regime, while we have a 
transition to Phantom Energy when u) < — 1. The particular case of uj = — 1 
corresponds to a cosmological constant. Nevertheless, neither Dark Energy nor 
Phantom Energy models appear to be satisfactory to explain the acceleration. 
A proposal to avoid Dark and Phantom energy comes form the so-called modi- 
fied gravity theories. In particular, one could consider the following replacement 
in the Einstein-Hilbert action[2 (k = 8irG) 



Abstract 




It is clear that other more complicated choices could be done in place of / (R) 4 . 
In particular, one could consider / (R, R^R^, R al 3 lS R al3li \ ■ ■ ■) or / {R, G) 
where G is the Gauss-Bonnet invariant or any combination of these quantities- 
One of the prerogatives of a / (R) theory is the explanation of the cosmological 
constant. Nevertheless, nothing forbids to consider a more general situation 
where a / (R) is combined with a cosmological constant A c , especially in the 
context of the Wheeler-DeWitt equation (WDW)[5j. For a / (R) = R, one gets 



H={2K) G ijkl TT»TT 
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where 3 R is the scalar curvature in three dimensions. The main reason to work 
with a WDW equation becomes more transparent if we formally re-write the 
WDW equation as [9] 
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where 



V = J d 3 xy/g 
is the volume of the hypersurface S and 

A s = {2k) G ijkl ir ij n kl - y/?R/ (2k) . 
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Eq. ([3]) represents the Sturm-Liouville problem associated with the cosmological 
constant. The related boundary conditions are dictated by the choice of the 
trial wavefunctionals which, in our case are of the Gaussian type. Different 
types of wavefunctionals correspond to different boundary conditions. We can 
gain more information if we consider gij — g^ + hij , where gij is the background 
metric and hij is a quantum fluctuation around the background. Thus Eq.f3]) 
can be expanded in terms of hij . Since the kinetic part of As is quadratic in the 
momenta, we only need to expand the three-scalar curvature J d 3 Xy/g 3 R up to 
the quadratic order. However, to proceed with the computation, we also need an 
orthogonal decomposition on the tangent space of 3- metric deformations [10l 1 1 1 j : 



h ij = -{ ( j + 2V-t i )g i j + {Lli) ij +hij. 



The operator L maps £j into symmetric tracefree tensors 
(LO ii = V i tj+V j Zi-~g ij (V-0, 



(6) 
(7) 



hjj is the traceless-transverse component of the perturbation (TT) , namely 

g*% = 0, V% = 



(8) 



1 For a recent riview on / (R), see Refs.[3||4], while a recent review on the problem of / (G) 
and / (R, G) can be found in Ref.[6][7]. 
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and h is the trace of hij. It is immediate to recognize that the trace element 
a = h — 2 (V • £) is gauge invariant. If we perform the same decomposition also 
on the momentum 7r lJ , up to second order Eq.Q becomes 



1 



d 3 a 



+ A| + AS 
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Concerning the measure appearing in Eq.Q, we have to note that the decom- 
position ([6]) induces the following transformation on the functional measure 
Vhij — > Vhj-jV^i'DaJi, where the Jacobian related to the gauge vector variable 
& is 

1. 



det Ag 



(10) 



This is nothing but the famous Faddev-Popov determinant. It becomes more 
transparent if £ a is further decomposed into a transverse part £j with V a £ J = 
and a longitudinal part {fH with £]j = V a ^>, then can be expressed by an 
upper triangular matrix for certain backgrounds (e.g. Schwarzschild in three 
dimensions). It is immediate to recognize that for an Einstein space in any 
dimension, cross terms vanish and J\ can be expressed by a block diagonal 
matrix. Since det AB = det A det B 
into 



the functional measure T>hij factorizes 
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with 



= Ag lJ — i? y acting on transverse vectors, which is the Faddeev- 
Popov determinant. In writing the functional measure Vhij, we have here ig- 
nored the appearance of a multiplicative anomaly [8 . Thus the inner product 
can be written as 



vhipevo^* [hj.] [e] m * [^] * [^] vi/ [(J ] 
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(12) 

Nevertheless, since there is no interaction between ghost fields and the other 
components of the perturbation at this level of approximation, the Jacobian 
appearing in the numerator and in the denominator simplify. The reason can 
be found in terms of connected and disconnected terms. The disconnected terms 
appear in the Faddeev-Popov determinant and these ones are not linked by the 
Gaussian integration. This means that disconnected terms in the numerator 
and the same ones appearing in the denominator cancel out. Therefore, Eq.Q 
factorizes into three pieces. The piece containing A^ is the contribution of 
the transverse-traceless tensors (TT): essentially is the graviton contribution 
representing true physical degrees of freedom. Regarding the vector term A|J, 
we observe that under the action of infinitesimal diffeomorphism generated by 
a vector field ej, the components of ([6|) transform as follows [TO] 



fr + 2V-£, 



hi 



(13) 
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The Killing vectors satisfying the condition Vi£j +Vj& = 0, do not change hij, 
and thus should be excluded from the gauge group. All other diffeomorphisms 
act on hij nontrivially. We need to fix the residual gauge freedom on the vector 
£j. The simplest choice is £, = O.This new gauge fixing produces the same 
Faddeev-Popov determinant connected to the Jacobian J\ and therefore will 
not contribute to the final value. We are left with 



(2) 
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+ 
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(ip<rW<r) 



(14) 

Note that in the expansion of J E cPx^fgR to second order, a coupling term 
between the TT component and scalar one remains. However, the Gaussian 
integration does not allow such a mixing which has to be introduced with an 
appropriate wave functional. Extracting the TT tensor contribution from Eq. ([3]) 
approximated to second order in perturbation of the spatial part of the metric 
into a background term §ij , and a perturbation hij , we get 



A^ = 



d 3 x^/§G 



ijkl 



(2k) K^ 1 - [x, x) ijU + -L (A L ) ° K± (x, x) lakl 



where 



(15) 
(16) 



is the modified Lichnerowicz operator and A x,is the Lichnerowicz operator de- 
fined by 

(A L h) t] = Ahij - 2Ri kjl h M + R lk h>* + R ]k h k ; A = -V Q V a . (17) 

Qtjki re p resen ts the inverse DeWitt metric and all indices run from one to three. 
Note that the term 

- mfrij + 3 Rhfj (is) 

disappears in four dimensions. The propagator K 1 - (x, x) iakl can be represented 
as 



K 1 - {~x,^y ) iak i =Y, 



^ )x (^)4l )x (y) 



(19) 



2A(r) 

where (~x) are the eigenfunctions of A L . r denotes a complete set of 

indices and A (r) are a set of variational parameters to be determined by the 
minimization of Eq. p5[) . The expectation value of A^ is easily obtained by 
inserting the form of the propagator into Eq. (|15[) and minimizing with respect 
to the variational function A (r) . Thus the total one loop energy density for TT 
tensors becomes 



A 
8^G 



2 ^ 



\f4 (t) + ^/w|(t) 



(20) 
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The above expression makes sense only for ojf (r) > 0, where u>i are the eigen- 
values of A L . Concerning the scalar contribution of Eq. (fT4]) , in Ref.[20] has 
been proved that the cosmological constant contribution is 



8^G 



u 2 {r) 



(21) 



where u) (r) is the eigenvalue of the scalar part of the perturbation. In the next 
section, we will explictly evaluate Eqs. (|20l 2ip for a specific background. 



2 One loop energy Regularization and Renor- 
malization for a / (R) = R theory 

If we consider a background of the form 

dv 2 

ds 2 = —N 2 (r) dt 2 + + r 2 (dd 2 + sin 2 9d<p 2 ) , (22) 

then, with the help of Regge and Wheeler representation [TJ], ^Ajh^ can be 
reduced to 

fi(x)=ul l f i (x) 1 = 1,2 , (23) 

where we have used reduced fields of the form fa (x) = F{ (x) / r and where we 
have defined two r-dependent effective masses m\ (r) and m 2 , (r) 

' m 2 (r) = ^(l- b -^)+4,b>(r)-4zb(r) 

(r = r(x)). (24) 

^ mI(r) = ^(l-^) + ^6'(r) + ^6(r) 

In order to use the WKB approximation, from Eq. (f23]) we can extract two r- 
dependent radial wave numbers 

k 2 (r,l,u,. hnl )=u; 2 nl - l -V±}l-m 2 (r) t = 1, 2 . (25) 

When b (r) = rt = 2MG, the effective masses can be approximated in the range 
where r € [rt, 5r t /2] with m\ (r) = — m\ (r) = ttIq (r). Such a restriction comes 
from the fact that the effective masses, in this range, represent short distance 
contribution. Indeed, we expect to receive large contribution from quantum 
fluctuations at short distances. It is now possible to explicitly evaluate Eq. ([2T)]) 
in terms of the effective mass. To further proceed we use the W.K.B. method 
used by 't Hooft in the brick wall problem[13] and we count the number of 
modes with frequency less than u)i, i = 1,2. This is given approximately by 

g(ui)= (™* u i (l,uj i )(2l + l)dl, (26) 
Jo 



dx 2 



1(1 + 1) 



(r) 
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where Vi (l,u)i), i = 1,2 is the number of nodes in the mode with (l,u>i), such 
that (r = r (x)) 

1 f + °° I 

Ui(l,Ui) = - dxyjk?(r,l,ui)- (27) 



Here it is understood that the integration with respect to x and ? ma x is taken 
over those values which satisfy fef (r, l,u>i) > 0, i = 1,2. With the help of 
Eqs. (|26l27p . Eq.® becomes 

A 1 JL» /-+ 00 



87rG 7T ^— ' / n gL> 



1=1 



This is the graviton contribution to the induced cosmological constant to one 
loop. The explicit evaluation of Eq. ([28|) gives 



^ i 2 r+°° 



= Pi + Pi = 



8nG 

1=1 



(29) 



where we have included an additional 47r coming from the angular integration. 
The use of the zeta function regularization method to compute the energy den- 
sities pi and p2 leads to 



P^ (e) 



raf (r) 
64tt 2 



e Vm-(r) 



1,2 , (30) 



where we have introduced the additional mass parameter p in order to restore 
the correct dimension for the regularized quantities. Such an arbitrary mass 
scale emerges unavoidably in any regularization scheme. The renormalization is 
performed via the absorption of the divergent part into the re-definition of the 
bare classical constant A, namely A — > Ao + A dlv . The remaining finite value 
for the cosmological constant reads 

^ = {pi[ii)+P2^))= P T Jf{H,T), (31) 

where pi (p) has the same form of pi (s) but without the divergence. The quan- 
tity in Eq. lpTj) depends on the arbitrary mass scale p. It is appropriate to use 
the renormalization group equation to eliminate such a dependence. To this 
aim, wc impose that [14] 

1 dA (p) d TT 
^G fX ^p- = V^ M - (32) 

Solving it we find that the renormalized constant Ao should be treated as a 
running one in the sense that it varies provided that the scale p is changing 

A (p,r) _ A (^ ,r) (r) p 

8irG ~ 8nG + 16tt 2 Mo' ^ 



G 



Substituting Eq.® into Eq.® we find 



A (l^o,r) 



8nG 32tt 2 



n4(r) 



In 



(34) 



If we go back and look at Eq. ((3|) , we note that what we have actually computed 
is the opposite of an effective potential (better an effective energy). Therefore, 
we expect to find physically acceptable solutions in proximity of the extrema. 
We find that Eq. ((3"4"l) has an extremum when 



1 = ml (r) Aq (mo,Q = Wq (r) = Mo 

e 4^2 " 8irG 64ir 2 4^ 2 e 2 ' 1 ' 

Actually Ao (/io, f) is a maximum, corresponding to a minimum of the effective 
energy. Note also that there exists another extremum when 

mo (r) = M = 0. (36) 

This solution corresponds to Minkowski space, producing no effect on the vac- 
uum. For this reason it will be discarded. On the other hand, the effect of 
the gravitational fluctuations is to shift the minimum of the effective energy 
away from the flat solution leading to an induced cosmological constant. If we 
apply the same procedure to the scalar part of the perturbation, we find that 
the only consistent solution is that A a — 0. Therefore, the whole contribution 
is due to the physical degrees of freedom: the gravitonjTBj. Plugging Eq. ([33|) 
into Eq. (|3"3"|) . we find 

Ao 0,f) _ A (no,r) , "i^(f),_ a* _™%(r)(, , A n 



ftrG SnG ' 16^ ln f = ^ [' + 41 Vj ' (3?) 

If we set fio = mp, where mp is the Planck mass, we can find that 

A ° ^ r) - when ji = exp (-^) Mo . (38) 



8nG n V 2. 

Nevertheless, jl is of the order of the Planck mass again, but unfortunately 
is a scale which is very far from the nowadays observations. However, it is 
interesting to note that this approach can be generalized by replacing the scalar 
curvature R with a generic function of R. Although a / (R) theory does not 
need a cosmological constant, rather it should explain it, we shall consider the 
following Lagrangian density describing a generic f(R) theory of gravity 

C = V=g{f{R)- 2A) , with f" ? 0, (39) 

where / (R) is an arbitrary smooth function of the scalar curvature and primes 
denote differentiation with respect to the scalar curvature. A cosmological term 
is added also in this case for the sake of generality, because in any case, Eq. (f39|) 
represents the most general lagrangian to examine. Obviously /" = corre- 
sponds to GR.[17j. The semi-classical procedure followed in this work relies 
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heavily on the formalism outlined in Refs.|20, 16 . The main effect of this re- 
placement is that at the scale fio , we have a shift of the old induced cosmological 
constant into 



Aq (mo,0 = 1 
8ttG y/h{Ej 



A (Mo,r) , 1 f j3 ^Rf'(R)-f(R) 



JJ2 



8ttG IQitGV va f {R) 



, (40) 



where V is the volume of the system. Note that when / (R) = R, consistently 
it is h(R) = l with 

We can always choose the form of / (R) in such a way Ao (/io, r). This implies 

AqCmqiI) _ _J L_ /■ ,3 X ^ Rf(R)-f(R) , 

A comment is in order. We have found that our calculation is in agreement with 
Ref. [TS] , where only the graviton contribution is fundamental. Note also the ab- 
sence of a Faddeev- Popov determinant. This is in agreement with Ref. [15] but 
also with Ref. |10j , where the Faddeev-Popov determinant appears when pertur- 
bations of the shift vectors are considered. The second comment regards our one 
loop computation which is deeply different form the one loop computation of 
Refs. |18U19| . where the analysis has been done expanding directly / (R). In our 
case, the expansion involves only the three dimensional scalar curvature. Note 
that with the metric (|22p and the effective masses (|24p. in principle, we can ex- 
amine every spherically symmetric metric. Note also the absence of boundary 
terms in the evaluation of the induced cosmological constant. 
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